syms a x n b t c sint lnz;
% 定义符号变量a、x、n、b、t、c、sint、lnz。

f = sym(a*x^n + b*t + c);
% 创建一个符号表达式f，表示a*x^n + b*t + c。

f1 = subs(f, 'a', sint);
% 将符号表达式f中的'a'替换为符号变量sint的值，并将结果存储在f1中。

f2 = subs(f1, 'b', lnz);
% 将符号表达式f1中的'b'替换为符号变量lnz的值，并将结果存储在f2中。

f3 = subs(f2, 'c', diff(exp(2*t)))
% 将符号表达式f2中的'c'替换为exp(2*t)的导数，并将结果存储在f3中。

f4 = subs(f, 'n', 3);
% 将符号表达式f中的'n'替换为数值3，并将结果存储在f4中。

f5 = subs(f4, 'c', pi)
% 将符号表达式f4中的'c'替换为数值pi，并将结果存储在f5中。

f6 = subs(f, 'c', 1:2:5)
% 将符号表达式f中的'c'替换为向量[1, 3, 5]，并将结果存储在f6中。

f7 = subs(f, 'c', [1, 3; 2, 4])
% 将符号表达式f中的'c'替换为矩阵[1, 3; 2, 4]，并将结果存储在f7中。
